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Abstract 

In a random unitary matrix model at large N, we study the properties of the ex- 
pectation value of the character of the unitary matrix in the rank k symmetric tensor 
representation. We address the problem of whether the standard semiclassical tech- 
nique for solving the model in the large N limit can be applied when the representation 
is very large, with k of order N. We find that the eigenvalues do indeed localize on an 
extremum of the effective potential; however, for finite but sufficiently large k/N, it is 
not possible to replace the discrete eigenvalue density with a continuous one. Nonethe- 
less, the expectation value of the character has a well-defined large ./V limit, and when 
the discreteness of the eigenvalues is properly accounted for, it shows an intriguing 
approximate periodicity as a function of k/N. 



1 Introduction and Results 



Unitary matrix models have many applications as a tool to study quantum systems, as well as 
an interesting mathematical structure [1]. They make an appearance in the two-dimensional 
lattice Yang-Mills theory that was solved in the large N limit by Gross and Witten [2] and 
they have a long history of applications to mean field models of lattice gauge theory [3-5] 
and lattice models with induced QCD [6-11]. They have interesting deformations which 
are related to integrable models [12-17] and they have recently been of interest as effective 
field theories for the low energy degrees of freedom in Yang-Mills theory defined on certain 
compact spaces [18, 19]. In the latter context, with the appropriate effective action, the 
matrix integral that we shall be interested in, shown in equation (1) below, computes the 
Polyakov loop expectation value in an effective theory corresponding to finite temperature, 
large N, four dimensional Yang-Mills theory where the three dimensional space is a sphere 
[20, 21]. 

In this Paper, we wish to point out an interesting behavior of the expectation values of 
the character of the unitary matrix for representations with center charge of order N, in 
the large N limit. We focus on the completely symmetric representations, denoted R = Sk, 
corresponding to a Young tableaux with a single row of k boxes. We will take k large, k ~ N, 
and typically in the region where k is somewhat larger than N. We strongly suspect that 
the behavior we find also occurs for other representations, but we shall not analyze them 
here. 1 The matrix integrals that we study have the generic form 

{TIrU) ~ J[dU]e-sP] ' W 

where U is an N x N unitary matrix, [dU] is the Haar measure for unitary matrices, S[U] 
is a class function of U (in that it obeys S[U] = S[VUV'] for any unitary matrix V) and R 
is an irreducible representation of SU(N). S[U] is of order N 2 in the sense that S[X] ~ N 2 
where X is the unit matrix, when N is large. 
For later reference, define T k by the relation 

(Tr Sk U) = e~ Nr * . (2) 

In the following, we shall consider a unitary matrix model with a generic action equipped 
with the properties that we have listed above, although in many instances, to be concrete, 
we will use the Gross- Witten model (3) since its eigenvalue density is known explicitly. In 
the Gross- Witten model [2], 

S[U] = [Tr U + Tr W] . (3) 
A 



1 Completely antisymmetric representations whose Young tableaux is a single column can also readily be 
analyzed and have some interesting phase structure [22], but due to the limit of N boxes on the length of a 
column in their Young Tableau, they do not display the behavior that we discuss here. 
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The matrix model in (1) is what is referred to as an "eigenvalue model" since it is possible 
to use the symmetry of the action to write the integral as an integral over diagonal matrices. 
Because it is a class function, S[U] is a function of the eigenvalues only, as is Tr R U. Therefore, 
the above matrix integral can be written in the standard way as an integral over the TV 
eigenvalues of the matrix U, {e 1 ^ 1 , . . . , e l ^ N }, 

IU<f>i |A(e^)| 2 e- 5 ^ ' { ) 

where |A(e l ^)| 2 = Yli<j l e "^ ~ e l ^ j \ 2 is the Vandermonde determinant. To evaluate this 
expression at large N, several assumptions are usually made: 

1. Saddle point approximation First there is the assumption that there is a good 
large N limit where the integral can be evaluated in a saddle point approximation. 
This includes the assumption that the effective potential <S e ff(0i) = S[(f>i] — In |A(e t,?;>i )| 2 
has isolated minima and the magnitude of the potential itself for generic values of fa 
is of order N 2 . Then, for small enough representations R, the expectation value in 
equation (4) can be computed by the saddle point approximation. We will denote the 
location of the minimum of S e g(<j)i) by (f>j and assume for simplicity that — n < 4>i < 
02 < • • • < (f>N < 7r. With U = diag(e^ 1 , . . . , e^ N ), we have simply 

\n(Tr R U) = InTr R U + C(l/iV 2 ) . (5) 



2. Probe approximation The saddle point approximation can still be applied to the 
integral at large N where the representation R is large, corresponding to a Young 
Tableau with k boxes and k is of order N, that is, ^ remains finite as iV — > oo [23- 
25, 21]. Moreover, to compute the leading term, which is of order N in ln(TrR[7), 
it is sufficient to assume that the position of the saddle point is not affected by the 
presence of an operator TtrII in the integrand. The upshot is that, for these large 
representations, (5) is replaced by 

\n(Tr R U) = InTr R U + 0(1) . (6) 

and has nontrivial content since ln(Tr#?7) is of order N . 

3. Continuum approximation An essential tool which is used to subsequently evaluate 
the matrix elements is the distribution function for the eigenvalues (pi. For large N, 
sums over eigenvalues are replaced by integrals over an eigenvalue density, 

= 4 £ W " #0 (7) 
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which is assumed to approach a piece-wise smooth function of in the limit where N 
is large. It is then employed to compute traces, for example, 

Tr U k = e ' kh ~ N I #p(</>K fe0 • (8) 
3 J 

In the following, we shall examine the reliability of the large N expansion for computing 
expectation values of the forms (Tr U k ) and (Trs k U) when k ~ N. We shall see that, for 
good reasons, the saddle point computation of (Tr U k ) fails once is large enough. On 
the other hand, the expectation value of the character, (Trs fc ?7), is much better behaved. 
For it, we shall find that points 1 and 2, which are essentially the saddle point and probe 
approximations, hold up to the largest values of -j| that we can study, and that (Tr Sk U) has 
a well-defined large N limit in the regime where is of order one. However, the assumption 
of a continuous eigenvalue distribution fails, leading to intriguing behaviour. 

To illustrate the latter point, let us assume that all aspects of the large N expansion 
are valid and use them to compute (Tr Sk U) in the large N limit. In that limit, we assume 
that the eigenvalues are classical, given by those values of <pi which minimize the effective 
potential, which we denote by 0j, and we can simply evaluate the expectation value by 
substitution: (Trs k U) = Tr,s fc £/ where U is the classical diagonal matrix. For concreteness, 
let us consider the strong coupling phase of the Gross- Witten model which has action S[U] 
given in (3). It is known to be solved by the eigenvalue density 

P(0) = ^(l + 2pcos(0)) , (9) 

where P = \ when A > 2. For the sake of this argument, we will assume that this eigenvalue 
density can be used to compute Tr U k . Using this density, traces are given by 




Tr U k = ^e lk ^ « I Np if \k\ = 1 , (10) 
j=i 1 if \k\ > 1 . 

The trace of a matrix in the irrep Sk can be written in terms of the eigenvalues as 

Tr 5fc t>= ex p(X>4>) , (11) 

jl<-<3k \o=l / 

which is equivalent, due to the Frobenius formula, to 

V=uE (flTruA , (12) 

' vts k \j=l J 
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where Sk is the symmetric group, the sum is over all its elements a, and where li, / 2 , . . ., l p 
are the lengths of the cycles in the decomposition into cycles of the element a. 
We can now combine equations (10) and (12) to obtain 

Tr Sk U = ^(Np) k ^((N/k)pe) h (13) 

or 

Tr Sk U « e~ fcln [^] . (14) 

Several things are interesting about this expression: as expected, since k ~ N, InTr s k U is 
of order N. For sufficiently large k/N, > ep, it is negative. This crossing from a scenario 
where Trs k U is an effectively infinite exponential of a positive quantity of order N to one 
which is exponentially small in iV was interpreted as a phase transition in references [20, 21]. 

To test the adequacy of this expression, we have performed several numerical computa- 
tions. In one, we treat U as a random variable, approximating the integrals in equation (4) 
by Gaussian integrals near the saddle point (the precise methodology is described in Section 
3). In another, we rely on both the saddle point approximation and the probe approxima- 
tion (but not the continuum approximation), treating U as a strictly classical variable but 
maintaining discreteness of the eigenvalue distribution. The results, at p — 0.45, are given 
in Figure 1. Three things are evident in this Figure: 

• There is a good match between the results of the saddle point computation plus con- 
tinuum distribution (14) and the numerics for smaller values of j^. This match ends 
abruptly at |f ~ 1.25. 

• For larger values of j|, the numerical computations agree with each other rather well, 
but disagree with (14). Agreement of the numerical computations suggests that the 
existence of a saddle point and use of the classical matrix U are still valid, apparently 
for the whole range of that we consider. However, the continuum approximation 
which led to the result (14) must fail for values at and above ~ 1.25. 

• Treating U as a random variable and treating it as a classical variable with the operator 
Trs k U in the probe approximation results in differences of sub-leading order, consistent 
with equation (6). 

Thus the saddle point approximation and the probe approximation appear to be valid for the 
computation of (Tr^C/), but the continuum approximation does not. Nonetheless, validity 
of the saddle point and probe approximations allow us to write, using equation (12), 

\ ' ^es fc \j=i / / ' aes k \j=i / 
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Figure 1: — 1\ = ln(Trs fe {7)/iV is plotted on the vertical axis as a function of k/N which is 
plotted on the horizontal axis, in the strongly coupled Gross- Witten model with p = 0.45 and 
N = 45. Crosses represent a computation where U is treated as a random variable, open 
circles represent a computation assuming the saddle point approximation and the probe 
approximation, i.e. treating U as a classical variable, but retaining discreteness of the 
eigenvalue distribution. The solid black line represents the result due to the continuum 
approximation, equation (14). 
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In the last equality, it is legitimate to substitute the classical discrete eigenvalues. 

We now draw the reader's attention to the strange feature of (Trs k U) as a function of k, 
namely the recurrence, or approximate periodicity, which is seen in Figure 1. The remainder 
of this Introduction will discuss the origin and interpretation of this feature, assuming the 
validity of the formula (15) above. 

Formula (15) implies that for the purpose of computing (Ti Sk U), we can treat U as a 
classical variable, frozen at the saddle point value U. To study the implications of this 
formula, we need to examine the properties of Tr U k for large k and large N, in the regime 
where k ~ N. We illustrate these properties in Figure 2. 

For p = 0.25, Figure 2(a), the results shown agree with the continuum approximation 
result (10) up to k/N ss 0.35. Above this point, Tr U k is no longer zero, as was predicted 
by equation (10), and is instead of order N (we should recall here that since U is a unitary 
matrix, |Tr U k \ < N, so above k/N w 0.35, Tr U k is quite large, reaching as much as half of 
its maximum value.). 

A reasoning for why this is so comes from the fact that, at any finite N, trace relations 
allow us to compute the traces of a single classical matrix Tr U k for k > N/2 from the first 
N/2 traces. It is therefore not possible for just the first trace to be arbitrary and non-zero, 
and the rest zero, as the continuum approximation would demand. 

Similarly, for general (but assumed smooth) eigenvalue distributions with a infinite num- 
ber of frequency components, Tr U k would have to decay to zero with k — > oo, which again 
does not appear compatible with trace relations. Thus, we might expect Tr U k to become 
noisy at some point in the vicinity of -j| ~ |, and this is indeed what we see for the case 
of p = 0.25 in Figure 2(a). Notice, in Figure 2(b), that at p = 0.45 which is close to the 
maximum value that p can reach in the strong coupling phase of the Gross- Witten model 
(before the phase transition to the weak coupling phase with a gapped distribution), the 
agreement with the continuum approximation is even weaker. 

More intuitively, consider the eigenvalues {0i, . . . , <Pn} and the spacing between consecu- 
tive eigenvalues — 0«. At large N, this spacing is approximately (Np^i))^ 1 where p(0) 
is the eigenvalue density. When summing over the phases e tk ^ 3 , we are probing the structure 
of the distribution of eigenvalues at wave-number k, and we cannot replace the discrete sum 
with an integral unless k- (spacing) <C 1, or k/N <C p. The continuum approximation to the 
density will thus hold for a finite k/N, as long as it is not too large, the cutoff being of order 
the inverse height of p(4>). The narrower the eigenvalue distribution, the larger the density 
p((f>) and therefore the continuum approximation will hold for a larger range of k/N. Only 
for the delta-function eigenvalue distribution will the continuum approximation hold for any 
k/N of order N°. 

Returning to the strongly coupled Gross- Witten model, the continuum approximation 
result (14), which is based on equation (10), implies that (Tr s k U) becomes very small once 
k/N is large enough. This is a generic feature of (Yi Sk U) which appears when Tr U k goes 
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Figure 2: Numerical simulation of the Gross- Witten model in the strong coupling phase with 
iV = 45 and (a) p = 0.25 and (b) p = 0.45. Tr U k is plotted on the vertical axis as a function 
of k/N on the horizontal axis. 
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Figure 3: A caricature of Figure 2. 

to zero fast enough as k — > oo. Define 7^ by (Tis k U) = exp(— JV^) in such a scenario. 
We can think of 7^ as as obtained in the continuum approximation and pictured with a 
solid line in Figure I. As we have just discussed, at any finite N, Tr U k cannot in fact go 
to zero for large k, due to restrictions placed on it by the trace relations. We will examine 
the consequences of this fact by assuming that Tr U k fails to be zero at precisely one large 
value of k, which we will denote with m, Tr U m — s ~ N. We will assume that m ~ N. 
This is motivated by Figure 2, where we see that Tr U k takes exceptionally large values at 
isolated points, the smallest of which is m = k ps 1.8 AT (for p = 0.45, Figure 2(b)). What 
we are doing is replacing the situation in Figure 2 with a caricature shown in Figure 3. (To 
be clear, 7*. could be computed by assuming that the only non-zero trace is Tr U ±l whereas 
Y k is computed assuming that the only nonzero traces are Tr C/ ±:L , Tr U ±m , though the 
argument is somewhat more general than that, as it relies on the large k behaviour of 7^, 
and not the details of Tr U k for small k.) 

Equation (12) now implies that, for n of order and k < m, 



Tr. 



5fc + 71 



u 



k\ 



(k + nm)\ 



(k + nm)\ 



a=0 



a\m a {k + (n — a)m)\ 



e - N lk + (n-a)r 



(16) 



Assuming that m is large enough that e N ^+nm <^ e N lk+{n _ l)m <^ <^ e Ny k ^ ^ e a = n 
term dominates the sum and we have 



k+nm n\m n 



1 



n\ \m/N 



(17) 



or 



T k+nm = T k + 0(\n(N)) . (18) 

We discover that the order N part of is periodic in k, with the period given by m. This 
simple calculation reproduces the behavior seen in Figure 1 quite well. In particular, we 



8 



see that the period of the recurrence in Figure 1 is about 1.8 which matches really well 
anomalously large value of Tr U k seen at k — 1.8N in Figure 2(b). 

Going back to Figure 1 we see that Tr U k is also large at k — 2m and k = 3m (and probably 
at higher multiples). How does that affect our computation? Let's say that Tr U 2m = s. 
Then, equation (17), for n — 2, becomes 

Tr — e 53 (h fe) 2 + it fey') w • (19) 

which gives T k+2m = Tfc + 0(1/ N) as before. Thus, the presence of further points where 
Tr U k is nonzero does not affect our conclusions as long as these appear at k equal to integer 
multiples of m. This turns out to be a generic feature of simple eigenvalue distributions, and 
we will discuss it further in Section 2. We will also examine a more complicated example 
with multiple periods in Section 4. 

Before moving on, let us briefly discuss the origin of the values of k for which Tr U k is 
anomalously large. They can occur when there is a region in the interval [— tt, it] where a 
large number of the classical eigenvalues are almost equally spaced. This would damp the 
destructive interference between the individual terms in the sum over j of the phases e lk( ^ j . 
It would occur for values of the wave- vector k which are 2n ■ integer/ (spacing). Eigenvalues 
are almost equally spaced when the derivative of the eigenvalue distribution is zero, that is, 
at an extremum of the distribution. We will return to this issue in the following Section and 
see a high degree of correlation between the period of recurrences and the maxima of the 
eigenvalue distribution. This would also explain why the anomalously large Tr U k seen in 
Figure 2 apparently occur for values of k which are integer multiples of a basic number. 

So far we have focused on the properties of Tr U k . What about (Tr U k )l Figure 4 
addresses this question. It shows (Tr U k ) in the Gross- Witten model with p = 0.25 computed 
by treating U as a random variable, correcting the classical limit by integrating over quadratic 
fluctuations around the saddle point (filled circles), as well the classical limit itself, Tr U k 
(from Figure 2(a)) for comparison. Curiously, the filled circles do agree with the large N 
result (10) over the entire range. We note that this integration is done for SU(N), rather 
than U (N) by imposing the constraint that the average value of the eigenvalues must be 
zero. If we use the unconstrained integral over U(N), there is a soft mode which changes 
the result considerably. We take the upshot of this discussion as evidence that the saddle 
point approximation cannot be trusted for the computation of (Tr U k ) when k is bigger than 
some fraction of N. The fluctuations are apparently at least as important as the classical 
eigenvalues in that regime. This is consistent with equations (5) and (6), since (Tr U k ) is 
bounded by N and therefore its logarithm cannot be large. The failure of the saddle point 
approximation arises from the variance in the eigenvalues fa being of the same order as the 
spacing between 0jS in the classical limit. Thus integrating over the eigenvalues 'washes' 
out the effect of a discrete spectrum and seems to restore the validity of the continuum 
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Figure 4: Numerical simulation of the Gross- Witten model in the strong coupling phase with 
p = 0.25 and iV = 45. (Tr U k ) is plotted on the vertical axis as a function of k/N on the 
horizontal axis. Open circles are obtained by finding 0j that minimize the effective action 
and evaluating Tr U k using those eigenvalues. This is identical to Figure 2(a). Filled circles 
include the leading correction to the large N result, arising from integrating the fluctuations 
about 4>i. The precise methodology is given in Section 3. 



10 



approximation for (Tr U k ). Another way to view this is that re-introducing some randomness 
to the matrix so that it is no longer strictly classical relaxes the trace relations. In fact, we 
see that this works rather well in this case, in that fluctuations allow the traces to agree 
with what one would obtain from the eigenvalue density in the large N limit over the entire 
range of that we explore. 

The curious (and central to our argument) fact is that this is not the case with (Trs k U). 
If we were to use (Tr Sk U) instead of Tr U k in equation (15), we would not obtain the correct 
answer. This can potentially be explained by a failure of factorisability (it is unlikely that 
Y\ k (Tr U ki ) = (J\ k Tr U ki ) if (Tr U k ) cannot be accurately computed in the saddle point 
approximation), but we do not explore this any further here. In contrast to (Tr U k ), (Tr Sk U) 
seems to be computable in the saddle point approximation (as long as we do not use the 
continuum approximation as well) and to have a well defined large N limit, which we will 
explore further in Section 2. 

The periodic behavior of equation (18) has an interesting implication for some physical 
applications of this matrix model. For example, when the character is the expectation value 
of the Polyakov loop in Yang-Mills theory on the sphere, T k is interpreted as the free energy 
of a heavy quark with center charge k and in the totally symmetric representation. If one 
thinks of this heavy quark as being composed of k partons, it appears that m of these partons 
can combine together to form a ('bound') state of low free energy. Thus, if k = nm + k 
(0 < k < m), the free energy of k partons receives contributions mainly from the k partons 
that are 'free', while the contribution of the n 'bound' states of m partons is negligible. We 
will argue that this periodicity always exists, since it is related to maxima of the eigenvalue 
distribution. 

The rest of the paper is organized as follows. In Section 2, we study further the character 
and its large N limit in the Gross- Witten model, both in strong and in weak coupling. 
We find that the character is indeed quasi-periodic and link the period to the maximum 
value taken by the eigenvalue density. In Section 3, we numerically compute the integrals in 
equation (4) to test approximations 1 and 2, and to see whether our results from Section 2 
are robust. Finally, in Section 4, we compute T k for more complicated eigenvalue densities 
and discuss generic behaviour. 

2 Character in the Gross- Witten model 

In this Section, we shall examine the difference between using a discrete and continuous 
eigenvalue distribution in the limit where N is large. We will assume that the large N 
limit localizes the matrix integral onto eigenvalues fa and that to a first approximation, 
these eigenvalues are not disturbed by the presence of the character in the integral. To aid 
numerical computation at large N, will use an approximation to the eigenvalues fa which we 



11 



1.4- 




Figure 5: Numerical simulation of the Gross- Witten model in both the strong and the weak 
coupling phases. <Jk on the vertical axis is plotted as a function of k/N on the horizontal axis 
for p=0.25, 0.40, 0.75, 0.90 (from the shortest recurrence period to the longest). Red dots 
represent numerical results at N=200 and the solid black line is the result in the continuum 
approximation. 
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obtain be re-discretizing the continuum distribution. We will show that this is in fact a good 
approximation in the next Section. Here, it will prove sufficient to demonstrate our point 
about the consequences of the failure of the continuum approximation. To be concrete, we 
will study the expectation value of a character in the Gross- Witten model [2], both below 
and above the third-order phase transition. 

At finite N, we approximate the actual positions of the eigenvalues, 0j, by approximate 
positions 4>i computed from the large- N density, via the formula which (at infinite N) defines 
the density, 

k-l 



f 



p(0) = for k = 1 . . . N . (20) 



If p(0) = p(—(f>), the | simply leads to a symmetric distribution of eigenvalues. 0jS differ 
from 0jS by corrections which are sub-leading in the large- N limit. We will demonstrate that 
they are a good enough approximation to (pi in the next Section. 

We will find it convenient to use the following generating function for the character: 

oo N 1 

e tk t ^ u = n t— ^- = ex p (- Tr ^ - tu )) ■ ( 2i ) 

k=0 j=l 

Then, a k defined by 

oo / n \ 

Y^t k e Nck = exp -iVE ln (l - te~*t) (22) 

k=0 \ j=l J 

is approximately equal to — T k . We will see in the next Section that the difference between 
(Tfe and — is of order 1/N. 

In the ungapped phase of the Gross- Witten model, the eigenvalue density is 

P(0) = ^(l + 2pcos0) (23) 

for < p < 1/2, while in the gapped phase, the eigenvalue density is 

rh 

COS 



M) = <2 l\ p) \J 2 ~ 2 P- sin2 (^/ 2 ) ( 24 ) 

for 1/2 < p < 1. Here, in both cases, we have parameterized the distributions by the 
expectation value p = ^(Tr U). With these explicit eigenvalue densities, at finite fixed 
N, we solve equation (20) (numerically) for 0j at different values of p. This allows us to 
write down the generating function (22) as a function of t, which we then expand (using 
Maple) in a Taylor series for small t to obtain a k as a function of k. The results, at N=200, 
are displayed in Figure 5. Since both axes of the Figure are scaled so that the plots are 
independent of iV (to leading order), finite iV results should approximate infinite iV results. 
For comparison, solid lines show the answer at infinite iV in the continuum approximation, 
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Figure 6: Numerical simulation of the Gross- Witten model for p = 0.55. on the vertical 
axis is plotted as a function of k/N on the horizontal axis for N=100 (green), N=200 (red) 
and N=400 (blue). The solid black line is the result in the continuum approximation. 

as computed in the introduction and in [20]. The approximate periodicity of as a function 
of k/N, predicted in the Introduction, is clearly visible in the Figure. 

One might worry that the recurrences shown in Figure 5 are simply a consequence of 
working at finite N and do not accurately represent the large iV limit of the theory. To show 
that this is not the case, Figure 6 contains the results of a computation at N = 100, 200 
and 400. Finite N results seem to converge to a well-defined answer with a stable recurrence 
period. 

We can understand the period of the recurrences as follows. In the Introduction, we 
argued that the recurrences are due to a single large trace, Tr U m at some m. Since Tr U m = 
J2j exp(zm0j), Tr U m is large if the spacing between consecutive <fij is an integer multiple of 
2n/m. Of course, the spacing is not constant, so what we want is that a large fraction of the 
eigenvalues be spaced at approximately 2n/m. The eigenvalue spacing, equal to (Np^cp))^ 1 
varies slowest where the derivative of p(<f>) is zero. Since there are more eigenvalues near the 
point where p(<p) attains its maximum than near the point where it attains its minimum, we 
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Figure 7: The recurrence period as a function of 2irp max . The crosses are data based on our 
numerical results and the solid line is our prediction that recurrence period = 2np max . 

conjecture that m is given by 2nN p maX) where p max is given by 



The recurrence period for as a function of k/N should then be simply 2irp max . This con- 
jecture is clearly supported by Figure 7. We will see in Section 4 that for more complicated 
eigenvalue densities with multiple maxima, there will be several competing recurrence peri- 
ods. Nonetheless, the periods will be related to the local maxima of the eigenvalue density 
in the way described above. 

Of course, the toy picture painted in the Introduction, where there Tr U m is large for only 
one isolated value of m is not accurate even in the case of the relatively simple eigenvalue 
distributions discussed in this Section. This was already illustrated in Figure 2. Here we 
complement that Figure with a computation in the weak coupling regime. Figure 8 shows 
Tr U k for p = 0.55 and N = 200. Several things are apparent in this plot: Tr U is large, while 
Tr U k is small at first, but grows with k. It also appears nearly random with a distribution 
whose mean is zero. Since Tr U k oscillates, the effect described in the Introduction from 
each individual Tr U k roughly cancels between different values of k. However, at k/N m 2.1 
something interesting happens: Tr U k is of order iV and negative for several values of k in 
a row. Thus, the effect for several values of k in a row can reinforce, which should lead to 
a recurrence. We see in Figure 6 that indeed, at p = 0.55, the recurrence period is about 2. 
Not surprisingly, a similar reinforcement will occur for k/N « 4.2. 




< p < 1/2 , 
1/2 <p< 1 . 



(25) 
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Figure 8: Tr U k in the Gross- Witten model for p = 0.55 and iV = 200. 

3 Matrix model beyond the saddle point approxima- 
tion 

In this Section we provide the methodology behind Figures 1 and 4 in the Introduction and 
further support the validity of our approximation in Section 2. While we have shown already 
that the recurrence pattern is stable in the large N limit, one could wonder whether it is also 
stable under small (of order l/N) corrections to the eigenvalue density. We have checked 
the impact of such order l/N changes as removing a single eigenvalue and adding a small 
high frequency component to the eigenvalue density, and these had no impact on the overall 
pattern beyond l/N corrections. 

Another possibility is that the since the eigenvalues are random variables with a variance 
of order l/N, the discreteness of the eigenvalue distribution might be washed out when 
the integrals in equation (4) are evaluated. To test for this possibility, in this Section we 
evaluate these integrals by approximating the potential near the saddle point (to the leading, 
quadratic order). While this computation is not exact, it should provide us with the next-to- 
leading order corrections to and allow us to see whether the discreteness of the eigenvalue 
density still matters when the eigenvalues are random variables. 
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The expectation value we are interested in can be written as 



(e 



-Trln(l- 



tU)j 



I Uf=i dx j e~ v ^ x ^ exp (- J^j ln(l - *e^)) 



where the effective potential for the eigenvalues is 



(26) 



Kff = -^ln 



i<j 



sin 



X r 



2pN ^ cos(xj) . 



(27) 



The first term in the potential, which comes from the Vandermonde determinant, causes 
eigenvalues to repel, while the second term attracts them towards Xj = 0. The second term 
comes from the Gross- Witten action (3) where A = 1/p for p < 1/2 and A = 4 — 4p for 
p > 1/2. This action is known[2] to produce the eigenvalue distributions discussed in the 
Section 2. 

We will denote by Xi the positions of the eigenvalues for which V c g attains its minimum. 
X{ are, up to corrections of order 1/N, the same as fa in equation (20). The minimum is 
unique up to a permutation of the eigenvalues x iy and satisfies ^Xi = 0. Given x iy we can 
expand V e s around its minimum 

V eS (x 1 ,...,x N ) « Vrff^i,...,^) (28) 

1 N d 2 

+ » X] q~ q~ V cS (x 1 ,...,x N )\ Xk=Xk (x a - x a ) (x b - x b ) . (29) 



2 ^ dx a dx b 

a, 6=1 



We also expand the inserted operator O = exp ^— ln(l — te lXj )*j around the point 

Xi Xi, 

N d 

0(x!,...,x N ) « 0(xi,...,x N ) + ^—0(xi,...,x N )\ Xk= x k (x a -x a ) (30) 



o=l 



+ 



1 N 



d 2 



2 ^ dx a dx b 

a, b=l 



0(xi, . . .,x N ) 



I x k — x k 



(x a -x a ) (x b -x b ) . (31) 



The second term on the first line does not contribute to the integrals, by symmetry. Evalu- 
ating the appropriate Gaussian integrals we obtain 



(32) 



where we have arranged the second-order derivatives into N x N matrices in the natural way. 
Xi, and the matrices of second-order derivatives can be obtained in Maple. The resulting 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 

k/N 

Figure 9: — Tj, in the Gross- Witten model as a function of k/N for p = 0.55. Open circles 
represent the saddle point result (the first term in equation (32)), while crosses include 1/N 
corrections (both terms in equation (32)); N=45. The black points show the result of an 
approximate saddle point computation with N=200 as discussed in Section 2, for comparison. 
The three plots differ by 1/N corrections. 

expectation value is a function of t, and can be expanded at small t as before. The first 
term in the above equation is simply the saddle point result we have obtained before (but 
with x instead of 0), while the second term represents 1/N corrections due to variance in the 
eigenvalues. As can be seen in Figure 9 (and in Figure 1), the general pattern we observed 
in the previous Section persists. The difference between the saddle point results at iV = 45 
(open circles) and at iV = 200 (points) is of order 1/N, as is the difference between the result 
of a Gaussian integral (crosses) and just the saddle point (open circles). This again justifies 
approximations 1 and 2 described in the Introduction, demonstrates that our methodology 
in Section 2 is sufficient and shows that the recurrences we have seen appear in the exact 
computation as well. 

To compute the results shown in Figure 4, we proceed in a similar fashion, except that we 
do not expand the operator, exp(ikxj) around Xj = Xj, as it is not necessary to do that to 
obtain the integrals. Also, we impose a constraint ^ X{ = in our integrals, corresponding 
to the SU(N) model, as the U(N) model has a soft mode which leads to (Tr U k ) ~ even 
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Figure 10: p((f>) given in equation (33) and the corresponding a k as a function of k/N, for 
N = 200. 

at small k (including k — 1). 



4 Further examples and discussion 

The results presented in Section 2 are quite generic if one focuses on eigenvalue distributions 
with a single maximum (either ungapped or one-cut gapped distributions). When multiple 
maxima are present, however, the behaviour can be much richer, as multiple maxima can lead 
to multiple recurrence periods. As an example, consider the following two-cut distribution, 
consisting of two rescaled semicircle pieces, one centered around = and one around = n 



P(0 = 



hi a /0 2 - wj 
h 2 y/(<f> ~ tt) 2 - w\ 



for |0| < w\ , 
for |0 — 7r| < w 2 , 



(33) 



where hi = 0.46283, hi = 0.20123, w\ = 0.097T and wi = 0.8n. The eigenvalue distribution 
together with the corresponding a k as a function of k/N are shown in Figure 10. Two 
different recurrence periods are clearly visible, each related to one of the two maxima by 
the formula (recurrence period) = 2np max . The eigenvalue distribution was carefully chosen 
to make the amplitude of the two series of recurrences approximately equal so that they 
would both be visible. Generically, even if the distribution has multiple maxima, one of the 
recurrence sequences dominates the others. 

Other phenomena are possible as well. For example, in Figure 11, with an ungapped 
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Figure 11: p(<f>) given in equation (34) and the corresponding a k as a function of /c/iV, for 
N = 200. 

distribution given by 



we see that the peaks can broaden and eventually merge as they reoccur. 

The recurrences as a function of k/N seem to be a fairly generic feature of T k at large 
N. This is interesting partly because for many theories of interest, such as Af = 4 SYM, the 
eigenvalue density is not known. Therefore, while we cannot compute in M = 4 SYM, 
we can make a conjecture that it will exhibit some kind of recurrences. Notice further that 
Tk seems also to generically experience intriguing first-order phase transitions (kinks) as a 
function of k/N. 

In the Introduction, we have already noted the qualitative difference between the be- 
haviour of (Tr Sk U) and (Tr U k ). Not only is the former exponentially larger than the latter, 
but it exhibits a very clean large N limit which can be obtained from the saddle point ap- 
proximation. In contrast, (Tr U k ) cannot be obtain this way. This suggests that if one is 
interested in studying the moments of the eigenvalue distribution, (Tr s k U) is a better object 
to compute than (Tr U k ). These two objects contain in principle the same information, 
but organized in a different way. In particular, through exponentiation, (Trs k U) magnifies 
certain features of the eigenvalue distribution p (such as the value of p ma x) making them 
easier to extract. 



p{<t>) 



1 



(-0.3 cos(0) + 0.15 cos(20) + 0.25 cos(30)) , 



(34) 
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Since (Tr Sk U) has a good large-N limit, it is also the natural candidate for computations 
using the AdS-CFT duality (the same is true of the expectation value of a character in 
the rank k totally anti-symmetric tensor representation, as was shown by similar methods in 
[22]). It would be most interesting to find a stringy counterpart of the recurrences conjectured 
in this Paper. 

Finally, the potentially complex recurrence pattern means that computing (Tr s k U) in 
finite temperature M = 4 SYM on a sphere (for example, through the AdS-CFT duality) 
might teach us about more than just the eigenvalue distribution. If the picture we presented 
in the Introduction, in which the recurrences are due to presence of bound states of partons, 
is correct, the recurrence pattern directly carries information about the dynamics of quarks 
in SYM. 
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